OH hyperfine ground state: from precision measurement to molecular qubits 
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We perform precision microwave spectroscopy — aided by Stark deceleration — to reveal the low 
magnetic held behavior of OH in its 2 n 3 /2 ro-vibronic ground state, identifying two held-insensitive 
hyperfine transitions suitable as qubits and determining a differential Lande p-factor of 1.267(5) x 
10~ 3 between opposite parity components of the A-doublet. The data are successfully modeled 
with an effective hyperfine Zeeman Hamiltonian, which we use to make a tenfold improvement of 
the magnetically sensitive, astrophysically important AF = ±1 satellite-line frequencies, yielding 
1 720 529 887(10) Hz and 1 612 230 825(15) Hz. 

PACS numbers: 33.20.Bx, 33.15.Pw, 33.55.Be, 39.10.+j 
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Scientific advances such as creating new forms of dipo- 
lar matter, realizing molecular scalable quantum com- 
putation, and uncovering physics beyond the Standard 
Model demand improved understanding of even the most 
familiar of polar molecules. Techniques such as buffer 
gas cooling Q , photoassociation , and Stark decelera- 
tion (3, |J] , are beginning to yield cold samples of ground 
state polar molecules, such as OH, that can be utilized 
for these endeavors. However, so far no technique simul- 
taneously achieves the large molecule number and ultra- 
cold, <100 /iK, temperature required. Solving this prob- 
lem would be exciting not only for the accompanying 
ability to study ultracold collisions and collective phe- 
nomena dominated by long-range dipole-dipole interac- 
tions 5], but also for the prospect of using the strong 
dipolar interaction to robustly implement quantum logic 
gates 0, 0|. This would obviate the need for the pre- 
cise, dynamical position control or the short-lived excited 
states that encumber other proposed quantum computa- 
tion systems 0. 

Stark deceleration currently provides relatively large 
numbers of polar molecules, but at temperatures limited 
to a few mK. Additional cooling will be required if we 
hope to reach densities necessary for exploring ultracold 
dipolar gas phenomena. Magnetic trapping is a promis- 
ing aid for cooling because it provides long term confine- 
ment of relatively hot particles while also admitting elec- 
tric (E) fields to the trapped molecules. This latter point 
is important since B- and E-fields can modify collisional 
processes in ways that enhance the efficiency of sympa- 
thetic cooling with co-trapped alkalis. Once cooled, E- 
fields can tune the dipole interaction to explore collective 
phenomena such as field- linked states yjg. To achieve 
these goals, we need to understand precisely the influ- 
ence of B-fields on collisional resonances in the cold and 
ultracold regimes, and experimental data are required to 
validate theory models of polar molecular states in low 
B-fields. Additionally, because fluctuating B-fields are 
a nuisance for long-term qubit coherence and precision 
measurements, accurate models are required to account 
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FIG. 1: (color online) (a) OH angular momentum perturba- 
tions. J, L, and S are the total (minus nuclear), orbital, and 
spin angular momentum operators, (b) Small B-field Zee- 
man shifts in 2 n 3/ / 2 - Long-dashed arrows depict main-line 
7r-transitions for SgF measurement while short-dashed show 
measured satellite lines, (c) Main-line 7r-transition frequency 
shifts, A/, versus B-field. Curves of identical dash pattern 
asymptote to same set of J-basis lines. Aa«1.7 GHz omitted 
from vertical scale. 



for and nullify these unwanted effects. The natural com- 
bination of precision microwave spectroscopy with the 
cold, monoenergetic packets from Stark decelerators || 
enables us to understand polar molecule hyperfine Zee- 
man behavior, as presented in this Letter. 

As a demonstration of this technique, we experimen- 
tally verify our refinement of previous OH Zeeman the- 
ory 0, Ho| to the hyperfine regime and use it to identify 
two candidate qubit transitions that are robust against 
B-ficld fluctuations. Moreover, we apply this model to 
the precision measurement of the magnetically sensitive 
OH ground state satellite-lines, completing the ten-fold 
resolution improvement in all four lines of the astro- 
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physically important J = 3/2 A-doublet These 
measurements highlight the ability of the Stark decel- 
erator not only to enhance our understanding of unex- 
plored regimes of molecular coupling, but also to con- 
tribute toward searches of non-Standard Model physics, 
e.g. the variation of fundamental constants such as fine- 
structure, a, and the electron-proton mass ratio, y, which 
may be measured by comparing Earth-bound OH with 
that found in OH megamasers 0, These dis- 

tant sources are spatially well-defined, and by combin- 
ing our recent measurements with astrophysical studies 
of comparable resolution, we will be able to constrain — 
with spatial dependence — fine-structure variation below 
1 ppm for Aa/a over 6.5 G yr 0|. We further note that 
the low-field behavior revealed in OH may also occur in 
other molecules and certainly in those with 2 H structure. 
Methods presented here are crucial if B-field effects are 
to be enhanced for molecular cooling schemes or miti- 
gated in molecular-based clocks, qubits [H0]> and preci- 
sion measurements of the electron electric dipole moment 
in molecules [lij . 

The hyperfine Zeeman structure of OH is subtly per- 
turbed by angular momentum couplings of higher order 
than the well-known, diagonal L • S spin-orbit interac- 
tion, and this study reveals their effect on OH tran- 
sition frequencies in low B-field. These perturbations 
cause the parity states, labeled / and e, of the 2 U 3 / 2 ro- 
vibronic ground state — which are otherwise identical — to 
possess slightly different hyperfine energy splittings and 
g-factors. Figure sketches the relevant energy lev- 
els and the higher-order angular momentum couplings 
between them. The two 2 U fine-structure manifolds are 
mixed with the first electronically excited state, 2 Si/2, 
via J • L and the off-diagonal L • S couplings. Mediated 
by these couplings, the spin-rotation induced splitting of 
the p-doubled 2 E 1 / 2 state causes the 2 n states to also 
split. The resulting A-doublet energy splitting, Aa, be- 
tween parity states is smaller in 2 n 3 / 2 due to the weaker, 
off-diagonal L • S coupling. The hyperfine parity differ- 
ence may be ascribed to the minute difference in coupling 
strength between 2 £1/2 and the energetically separated 
parity states. Likewise, the larger /-state g-factor is due 
to the stronger J • S coupling between the energetically 
closer 2 n /-parity states |l5j. 

Experimentally, we focus on the electric dipole allowed, 
7r (Amp—0) and a (Amjf=±l) transitions between the 
A-doublet parity states, as shown by arrows in Fig.^b). 
We write these states in the \F,itif,p) basis and distin- 
guish 7r-transitions from states by substituting Ap for 
the parity value, p. The hyperfine energies of the parity 
states are written as F*' e = f ± £/2, where the average 
splitting, f, and difference, £, are approximately 54.1 and 
1.96 MHz, respectively |l6j. Similarly, the g-factors are 
written gf< e = g ± Sg/2, where Sg is 0.3% of g 

This work distinguishes itself from earlier studies [ToL 
ll^ l in which relatively poor spectroscopic resolution 
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FIG. 2: (color online) Data and theory fits for energy dif- 
ferences between main-line |2,mF,Ap) 7r-transitions versus 
B-field. Inset shows summary of Sg data. Red (large) hatch 
is previous 5gj data [To||. and blue (small) hatch is SgF- 

of the ground state precluded low- field measurements. 
While these subtle differences can be neglected in pre- 
vious experiments wherein huge magnetic fields cause 
the Zeeman energies to be much larger than the hy- 
perfine splittings, we demonstrate that — in this new era 
of high-resolution molecular spectroscopy — these effects 
cannot be ignored when low- fields are present. In other 
words, in fields of 1 T, where J is a good quantum num- 
ber, A-doublet transition frequency differences vary as 
expected — linearly with applied B-field |l0| . However, in 
fields <30 G, where F — J ± 1/2 is a good quantum 
number, we do not find a linear variation of transition 
frequency with B-field. We instead observe a linear fre- 
quency shift out to ~1 G, and thereafter detect shifts 
that are significantly varied in curvature versus B-field, 
as shown in the data of Fig. [21 for B < 8 G. 

Qualitatively, this effect may be understood as a differ- 
ential onset of the Paschen-Back effect in the two parity 
states, which we explicate by first considering the hypo- 
thetical case in which no differences exist between either 
the g-factors (Sg = 0) or the hyperfine splittings (£ = 0). 
As such, the parity states are identical with respect to 
an applied B-field, and one would expect no frequency 
shifts in the hyperfine Zeeman transitions of, e.g., the 7r- 
transitions between main-lines (AF = 0) versus B-field 
until B> IT. This holds in spite of quadratic level shifts 
induced by mixing between F = 2 and F — 1 at B > 30 
G, as described by the Paschen-Back effect, since this on- 
set occurs at exactly the same B-field for p = /, e. If we 
now account for Sg ^ 0, then a linear variation of the hy- 
perfine Zeeman transitions appears and is proportional to 
mpSg at low fields. Finally adding the hyperfine pertur- 
bation induces an interesting dynamic in the frequency 
differences. The onset of the Paschen-Back effect now 
occurs at slightly different B-fields for the parity states, 
causing the otherwise linear differential frequency shifts 
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to become quadratic at small B-ficlds for transitions in- 
volving | mj? | = 1. 

To investigate this more rigorously, we employ an effec- 
tive Hamiltonian 0] to systematically distill the full — 
though intractable — system Hamiltonian into one that 
focuses solely on the state of interest, 2 n 3 / 2 - The re- 
sulting Zeeman Hamilitonian is expanded into the hy- 
pcrfine basis, and the hyperfine energy and g- factor dif- 
ferences are incorporated parametric ally , a point of de- 
parture from previous treatments [3, llfij . The resulting 
hyperfine Hamiltonian is applicable to higher-rotational 
states (J > 3/2) as well as to the 2 H 1 / 2 manifold (is| : 
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F = J+l/2 F=J-l/2 
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where (3^' e = ^Bg^' e . Equation represents two par- 
ity Hamiltonians for each value of mp, and upon di- 
agonalization yields 16 energy eigenvalues for the four 
J = 3/2 hyperfine states. The off-diagonal terms vanish 
for \mp\ = 2 as expected for states which do not couple 
to F = 1, and the diagonal Zeeman terms vanish for the 
linearly insensitive mp = states. Eigenvalue differences 
correspond to transition energies as a function of B-field. 

The main-line 7r-transitions used to investigate this 
model are depicted in Fig. QJb) and plotted in Fig. QJc) 
as a function of B-field using the eigenvalues from Eq. ^ 
All lines — except those with \rap\ = 2 — significantly de- 
part from linearity before they asymptote to their J-basis 
behavior. Of note is the |2,l,Ap) line which exhibits 
a negative frequency shift from its zero field value over 
a region spanning nearly 500 G. Figure [21 displays the 
F = 2 7r-transitions measured in the low-field, B < 10 
G regime. Evidently, the |2, l,Ap) line crosses zero fre- 
quency shift twice before pairing with its J-basis partner, 
1 2, 2, Ap), at high B-field. Interplay between two small 
quantities — £ and Sg — governs this behavior, which may 
be seen from an expansion of the eigenvalue energy dif- 
ferences to third-order in B-field, valid up to ^30 G: 
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A practical consequence of this interplay is the insen- 
sitivity to magnetic fluctuations of the |2, 1, Ap) and 
|1,— l,Ap) lines, indicated in Figs. |21 and [He), at B = 
3.01(2) G and 7.53(4) G, respectively. With experimen- 
tally achievable 1 mG B-field stability, we estimate fre- 
quency fluctuations to be <5 mHz, which provides a tran- 
sition stability superior to that of the rap — satellite- 
line components in zero field. In addition, each is a tran- 
sition between states of equal, non-zero magnetic mo- 
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FIG. 3: (a) Rabi spectroscoy of the |2, tuf, Ap) main line at 
B = 2.96 G. (b) Fit to |2, 0, /) -» |1,0, e) satellite line. Each 
data point is the weighted average of several measurements 
taken over many days. 



ment, allowing stable magnetic trapping without heat- 
ing. These states have large, equal-but-opposite electric 
dipoles = 1.67 D) which are easily aligned with low 
E-fields (~1 kV/cm), and as such are enticing candidates 
for encoding B-field stable molecular qubits [6j. 

To experimentally validate our model, we utilize the 
previously discussed Stark decelerator and microwave 
spectroscopy |8j to obtain high spectral resolutions in the 
presence of <10 G magnetic fields. In brief, the Stark 
decelerator provides OH primarily in the strongest, weak 
electric (E) field-seeking states, \F,mp,p) — |2, ^0, /). 
Typical operation produces a 10 6 cm -3 dense packet with 
mean velocity controllable from 410 m/s down to nearly 
rest and a minimum temperature of 5 mK |4j . We choose 
to run the decelerator at 200 m/s as a compromise be- 
tween slowing efficiency, microwave power, and interro- 
gation time We perform Rabi spectroscopy by in- 
terrogating the OH packet with one or more microwave 
pulses referenced to the Cs standard |8(. Laser- induced 
fluorescence provides state detection by exciting the / 
manifold — equally driving the two hyperfine states — with 
282 nm light and registering a fraction of the subsequent 
313 nm emission with a photomultiplier tube. Main and 
satellite transitions originating from |2, mp, /) are driven 
with a single microwave pulse and detected by registering 
a reduction in photon counts (see Fig^b)). 

The 10 cm long, TMoio microwave cavity is tuned near 
the |2, mp,Ap) main line. A solenoid is wound around 
the cavity and encased in a p.- metal magnetic shield. The 
latter reduces the ambient field to <6 mG, while the for- 
mer can apply up to a 10 G field along the resonant 
TMoio axis, which itself is aligned along the OH beam 
path. At the \2,mp 1 Ap) frequency, the cavity E- and 
B-ficlds are nearly collinear — determined by the absence 
of observed er-transitions — and the E-field magnitude is 
constant over 80% of the cavity. In contrast, the satellite- 
line frequencies are detuned far enough from cavity reso- 
nance to distort the associated intracavity E-field magni- 
tude and direction, and we observe both tt- (E || B) and 
a- (E _L B) transitions on these lines. 

We calibrate the intracavity B-field by tracking the 
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\2,rriF,f) — > |1, m' F , e) satellite splitting as a function of 
B-field in the <1 G regime where the Zeeman theory is 
well-understood. Concurrent observation of the n- and 
CT-lines provides cross-checks on the B-field calibration. 
Their relative peak heights and frequency shifts versus 
cavity position also provide a tool for mapping E-field 
direction and power inhomogeneities. When combined 
with the \2,mp, Ap) data, this mapping provides a mea- 
surement of the B-field inhomogeneity as well. We per- 
form B-field calibration and corresponding Zeeman shift 
measurements over several portions of the cavity. This 
reduces systematic errors, and we conclude that our B- 
field calibration and inhomogeneity mapping are accurate 
to within less than 1%. Other systematic shifts, such 
as from collisions, stray electric field, dopplcr shift, and 
black body radiation, are collectively <3 Hz 

Figure OJa) shows typical main-line data. Only four 
peaks are visible since for 7r-transitions mp — m' F = 
is selection- rule forbidden when AF = 0. The data are 
fit with four Rabi peaks of the form sin 2 (^^), where 
lu' r = y'wl + S 2 , and uir and 5 are the Rabi frequency 
and detuning, respectively. Each line's detuning is a free 
parameter. Figure [21 summarizes results taken at differ- 
ent B-fields. The larger error bars are from traces taken 
at microwave interrogation times as short as t — 70 fis, 
which have larger transit-limited resolutions than those 
taken at 400 fis. As expected, the tuf — m' F = ±2 lines 
are linear and we detect curvature in the mp = m' F = ±f 
lines at fields greater than ~1 G. The theory and data 
show excellent agreement without any free parameters. 
The only parameter in Eq. ^ with uncertainty greater 
than ^10 -5 is 5g, which was measured by Radford in 
the J-basis at B = 0.6-0.9 T to be Sgj = 1.29(2) x 10~ 3 . 
Coupling to higher rotational states at large fields may 
inflate the gj difference over that measured in the _F-basis 
regime. To test this, we perform a single-parameter fit 
to our low-field data to obtain SgF. The resulting fits 
are shown in Fig. |5J which are indistinguishable from the 
theory curves. The inset plots the measured SgF for each 
of the four mp states, along with their weighted mean, 
Sg F = 1.267(5) x 10 -3 . This value is la consistent with 
Radford's gj, but higher precision measurements in both 
regimes may distinguish a difference. 

As a further application of this hyperfine Zeeman 
theory for OH, we measure the magnetically sensitive 
satellite-line transition frequencies. Specifically, we trace 
the ^ 2 i = [2,0,/) -► |l,0,e) and v 12 = |2,0,e) |1,0, /) 
frequencies versus B-field. These Zeeman transitions are 
first-order insensitive to B-field, which suppresses spu- 
rious frequency shift contributions from calibration and 
ambient field uncertainties. Figure GJb) shows the ap- 
proximately quadratic variance of V21 versus applied B- 
field. Nearly overlapping tt- and er-transitions preclude 
an accurate measurement at zero applied field, and in- 
stead we must magnetically split the lines. The theory 
fit is used to project the data down to zero applied B- 



field where the ambient field is the only magnetic per- 
turbation. The uncertainty due to the ambient field on 
the line frequencies is estimated to be <0.6 Hz. The fit 
value for the B-field calibration is consistent with that ob- 
tained earlier, albeit with a larger uncertainty. We choose 
to use the fit-value's error to obtain a conservative esti- 
mate of the transition measurement's uncertainty. The 
fit has a reduced \ 2 — 1-2, and projecting the curve 
to zero applied field results in a bare frequency of 1/21 — 
1720529887(10) Hz. For the v 12 line, we use the distorted 
E-field at the cavity entrance to drive |2, ±1, /)— >|2, 0, e) 
cr-transitions. A second microwave pulse then drives the 
7r-transition to |1,0,/). Using the previous fit's calibra- 
tions, we measure v\2 = 1612 230825(15) Hz as the tran- 
sition's bare frequency. These values provide a tenfold 
improvement over the previous measurement accuracy of 
100-200 Hz 0. 

Variation of a may be constrained by taking the sum 
and difference of either the main- or satellite-lines. The 
difference is sensitive to a while the sum measures the 
red-shift systematic. The satellite-lines offer a 55-fold in- 
crease in resolution over the main-lines due to the 2f 
versus £ separation [ll|. Now that the complete A- 
doublet is measured with high-precision, one can simulta- 
neously constrain variation of a and y [l2j . All four lines 
must be detected from the same source, and the closure 
criterion, i.e., the zero difference of the averages of the 
main- and satellite-line frequencies, provides a critical 
systematics check. We obtain an Earth-bound closure 
of 44±21 Hz. Imminent astrophysical measurements of 
the satellite-lines below the 100 Hz level 0] highlight 
the prescience of our study for observing the variation of 
fundamental constants. 
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